DISPERSIVE PROPERTIES FOR DISCRETE SCHRODINGER 

EQUATIONS 

LIVIU I. IGNAT, DIANA STAN 

Abstract. In this paper we prove dispersive estimates for the system formed by two 
coupled discrete Schrodinger equations. We obtain estimates for the resolvent of the 
discrete operator and prove that it satisfies the limiting absorption principle. The decay 
of the solutions is proved by using classical and some new results on oscillatory integrals. 
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1. Introduction 
Let us consider the linear Schrodinger equation (LSE): 



< 

-h ; , , f iu t + u xx = 0, x 6l, t^O, 

1 > \ u(0,x) =<p(x), xeR. 

Linear equation (11.11) is solved by u(t,x) = S(t)ip, where S(t) = e ltA is the free 

Schrodinger operator. The linear semigroup has two important properties. First, the 

conservation of the L 2 -norm: 
> 

(I- 2 ) \\S(t)<p\\ L 2(R) = IMU 2 (M) 

m ■ and a dispersive estimate of the form: 

^: i 

(1-3) \(S(t)<p)(x)\ < (47r|t|)1/2 lMUi (R) , x G R, t ± 0. 

The space-time estimate 

(1-4) \\S(-)(p\\ L e(R,Le(m:)) < C'IMIz^r), 

due to Strichartz p2], is deeper. It guarantees that the solutions of system ( 11. ip decay 
as t becomes large and that they gain some spatial integrability. Inequality (11. 4p was 
generalized by Ginibre and Velo [3j. They proved the mixed space-time estimates, well 
known as Strichartz estimates: 

(1-5) \\S(-)(p\\ L <i(R,Lr(R)) <C(q,r)\\(f\\ L 2 m 

for the so-called admissible pairs (q,r): 

(1.6) _ = _(___), 2<g,r<oo. 
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Similar results can be stated in any space dimension but it is beyond the scope of this 
article. These estimates have been successfully applied to obtain well-posedness results for 
the nonlinear Schrodinger equation (see [2], [Ej and the reference therein). 
Let us now consider the following system of difference equations 

f iu t + A d u = 0, j ez,t^o, 

(L7) 

I w(0) = <p, 

where A^ is the discrete laplacian defined by 

(A d u)(j) = u j+l - 2uj + Uj-x, j E Z. 

Concerning the long time behavior of the solutions of system (11.71) in [IT] the authors have 
proved that a decay property similar to the one obtained for the continuous Schrodinger 
equation holds: 

(1.8) \\u{t)\\ la o {Ij) <C{\t\ + l)- l ^y\\ Hlj) , Vt^O. 

The proof of (jl.8p consists in writing the solution u of (11.7J) as the convolution between a 
kernel K t and the initial data (p and then estimate K t by using Van der Corput's lemma. 
For the linear semigroup exp(itA d ), Strichartz like estimates similar to those in (11.51) have 
been obtained in [TTJ for a larger class of pairs (q, r): 

(1-9) -<i(J--), 2<q,r<oo. 

q S y 2 r 

We also mention [5] and [6] where the authors consider a similar equation on KL by 
replacing A^ by Ad/h 2 and analyze the same properties in the context of numerical ap- 
proximations of the linear and nonlinear Schrodinger equation. 

A more thorough analysis has been done in [9] and [10] where the authors analyze the 
decay properties of the solutions of equation iu t + Au = where A = Ad — V , with 
V a real- valued potential. In these papers / 1 (Z) — /°°(Z) and l 2 _ a {1) — Z^(Z) estimates 
for exp(itA)P atC (A) have been obtained where P a>c (A) is the spectral projection to the 
absolutely continuous spectrum of A and /± CT (Z) are weighted Z 2 (Z)-spaces. 

In what concerns the Schodinger equation with variable coefficients we mention the 
results of Banica [JJ. Consider a partition of the real axis as follows: — oo = Xq < X\ < 
• • ■ < x n+ i = oo and a step function a(x) = bj 2 for x E (xi,x i+ i), where bi are positive 
numbers. The solution u of the Schrodinger equation 

iut(t, x) + (cr(x)u x ) x (t, x) = 0, for x E M, t ^ 0, 

u(0,x) = uo(x), ieR, 

satisfies the dispersion inequality 

||m(*)IU»(r) < C\t\- l/2 \\u \\ L i m , t ^ 0, 

where constant C depends on n and on sequence {&j}™ =0 - We recall that in [1] the above 
result was used in the analysis of the long time behavior of the solutions of the linear 
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Schodinger equation on regular trees. In the case of discrete equations the corresponding 
model is given by 

{LW) I U(Q) = <p, 

where the infinite matrix A is symmetric with a finite number of diagonals nonidentically 
vanishing. Once a result similar to [1J will be obtained for discrete Schrodinger equations 
with non-constant coefficients we can apply it to obtain dispersive estimates for discrete 
Schrodinger equations on trees. But as far as we know the study of the decay properties 
of solutions of system (11.101) in terms of the properties of A is a difficult task and we try 
to give here a partial answer to this problem. In the case when A is a diagonal matrix 
these properties are easily obtained by using the Fourier transform and classical estimates 
for oscillatory integrals. 

The main goal of this article is to analyze a simplified model which consists in coupling 
two DSE by Kirchhoff's type condition: 

iut(t,j) + b± 2 (A d u)(t,j) = j<-l,t^0, 



;i.ii) 



iv t {t,j) + bz 2 {A d v)(t,j) = 


J>1, 


u(t,0) =v(t,0), 


t^O, 


b^ 2 (u(t,-l)-u(t,0)) = b^(v(t,0)- 


-v{t,l)), t^O, 


u(0J) = ip(j), 


J<-1 


v(0,j) = <p(j)> 


J>1- 



In the above system u(t, 0) and v(t, 0) have been artificially introduced to couple the two 
equations on positive and negative integers. The third condition in the above system 
requires continuity along the interface j = and the fourth one can be interpreted as the 
continuity of the flux along the interface. 

The main result of this paper is given in the following theorem. 

Theorem 1.1. For any if G Z 2 (Z\ {0}) there exists a unique solution (u,v) G C(R, Z 2 (Z\ 
{0})) of system (II. lip . Moreover, there exists a positive constant C(&i,&2) such that 

(1.12) ||(u,v)(<)||i-(z\{o}) < C(6i,&2)(|*| + l)- 1/3 ||^||^(z\ { o}), Vt G R, 

holds for all cp e I 1 {Z\{0}). 

Using the well-known results of Keel and Tao [7] we obtain the following Strichartz-like 
estimates for the solutions of system (11.111) . 

Theorem 1.2. For any if G / 2 (Z\ {0}) the solution (u,v) of system (II. lip satisfies 

||(«,'y)IU?(K,r(z\{0})) < C(q,r)\\<p\\pp,\{o}) 
for all pairs (q,r) satisfying (11.91) . 
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The paper is organized as follows: In section [2] we present some discrete models, in 
particular system ( II. lip in the case 61 = 62 an d show how it is related with problem (jl.7p . 
In addition, a system with a dynamic coupling along the interface is presented. In section 
[3] we present some classical results on oscillatory integrals and make some improvements 
that we will need in the proof of Theorem ll.il In section H] we obtain an explicit formula for 
the resolvent associated with system (11.111) . We prove a limiting absorption principle and 
we give the proof of the main result of this paper. Finally we present some open problems. 

2. Some discrete models 

In this section in order to emphasize the main differences and difficulties with respect to 
the continuous case when we deal with discrete systems we will consider two models. In 
the first case we consider system (11.1 ip with the two coefficients in the front of the discrete 
laplacian equal. In the following we denote Z* = Z \ {0}. 

Theorem 2.1. Let us assume that &i = b 2 . For any ip G / 2 (Z*) there exists a unique 
solution u G C(R, Z 2 (Z*)) of system (ll.lip . Moreover there exists a positive constant C(b\) 
such that 

(2.1) IK*)lli«(z.) < c(6i)(|t| + rr^iMi^*), v*6i, 

holds for all (pel 1 {%*). 

In the particular case considered here we can reduce the proof of the dispersive estimate 
(12. ip to the analysis of two problems: one with Dirichlet's boundary condition and another 
one with a discrete Neumann's boundary condition. 

Before starting the proof of Theorem 12 .11 let us recall that in the case of system (11.71) its 
solution is given by u(t) = K t * (p where * is the standard convolution on Z and 

g-Ktrfn'CDgij^ t G M, J G Z. 

In [TT] a simple argument based on Van der Corput's lemma has been used to show that 
for any real number t the following holds: 

(2.2) ^(j^Cdtl + l)- 1 / 3 , Vj'GZ. 

Proof of Theorem \2.1[ The existence of the solutions is immediate since operator A defined 
in (12. 7p is bounded in / 2 (Z*). We prove now the decay property (12. ip . Let us restrict for 
simplicity to the case b\ — 62 — 1- 

For [u, v) solution of system (11.111) let us set 

SU) = vU) + «(-J) | D{j) = m^z2l j>o. 

Observe that u and v can be recovered from S and D as follows 

(u,v) = ((S-D)(-),S + D). 
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Writing the equations satisfied by u and v we obtain that D and S solve two discrete 
Schrodinger equations on Z + = {j G Z, j > 1} with Dirichlet, respectively Neumann 
boundary conditions: 

( iD t (t,j) + (A d D)(t,j) = j>l,t?0, 

(2.3) I D(t,0) = 0, t^O, 

[ D(0,j) = *®-p- j \ j>l, 

and 

( iS t (t,j) + (A d 3)(t,j) = Q j>l,t?0, 

(2.4) J S(t,0) = S(t,l), MO, 

Making an odd extension of the function D and using the representation formula for the 
solutions of ( 11.71) we obtain that the solution of the Dirichlet problem (12.31) satisfies 

(2.5) D(t,j)=^2{K t {j-k)-K t {j + k))D(0,k), t^0,j>l. 

k>\ 

A similar even extension of function S permits us to obtain the explicit formula for the 
solution of the Neumann problem (12.41) 



(2.6) 



Sit, 3) 



fc>i 



(K t (k-j) + K t (k + j-l))S(0,k), t^O, j>l. 



Using the decay of the kernel K t given by (j2.2p we obtain that S[t) and D{t) decay as 
(|t| + l) -1 / 3 and then the same property holds for u and v. This finishes the proof of this 
particular case. □ 

Observe that our proof has taken into account the particular structure of the equations. 
When the coefficients b\ and 62 are not equal we cannot write an equation verified by 
functions D or S. 

We now write system (jl.lip in matrix formulation. Using the coupling conditions at 
j = system (jl.lip can be written in the following equivalent form 

iU t + AU = 0, 
17(0) = <p, 

{u,v) T , u = (u0'))j<-i> v = ( v i)j>i and 

/ 

K 2 



where U 



(2.7) 



A 



bf 







V 



-267/ 

K 2 








-K 



W+% 














&f+&| 



b~ 2 2 



6 2 " 2 







bf 

-26 2 " 2 







T 2 



\ 






... y 
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1 the operator A can be decomposed as follows 
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In the particular case b\ = b 2 



A = A,+B 



However, we do not know how to use the dispersive properties of exp(ztArf) and the par- 
ticular structure of B in order to obtain the decay of the new semigroup exp(it(Ad + B)). 
Another model of interest is the following one inspired in the numerical approximations 
of LSE. Set 

/ x / ^r 2 , x < °> 

a(x) = < ,i 2 r, 

Using the following discrete derivative operator 

(du)(x) = u(x + — ) — u(x — -) 
we can introduce the second order discrete operator 

d(adu)(j) = a(j + -)u(j + 1) - (a(j + -) + a(j - -))u{j) + a(j - -)u(j - 1), j G Z. 



In this case we have to analyze the following system 
f tu t {t,j) + b^{A d u){t,j) = 0, 
iu t (t,j) + b 2 - 2 (A d u)(t,j) = 0, 

iu t {t, 0) + b^ 2 u(t, -1) - (b^ 2 + b 2 2 )u(t, 0) + b 2 l u{t } 1) 
k u{0,j) = tp(j), 



(2.8) 



J < -i, t ^ 0, 
J>1, Mo, 
M0, 

jeZ. 



In matrix formulation it reads iU t 
given by the following one 

/ 



AU = where U = (wfj))^, and the operator A is 



(2.9) 



.4 



-2foi 

o 





-2 





K 2 
iK 2 ' 








b- 2 2 





b- 2 2 ) 



b 2 2 

-2K 2 









b- 2 



o \ 








bf 





\ 

Observe that in the case b\ = b 2 the results of [11] give us the decay of the solutions. 

Regarding the long time behavior of the solutions of system (12.81) we have the following 
result. 



/ 



Theorem 2.2. For any if G / 2 (Z) there exists a unique solution u G C(M, / 2 (Z)) of system 
(12. 8p . Moreover, there exists a positive constant C(bi,b 2 ) such that 

i»(Z) <C(6i,6 2 )(|£| + ir 1/3 |MI^)> 



u 



Vte 
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holds for all tp G / 1 (Z). 

The proof of this result is similar to the one of Theorem 11.11 and we will only sketch it 
at the end of Section HI 

3. Oscillatory integrals 

In this section we present some classical tools for oscillatory integrals and we give an 
improvement of Van der Corput's Lemma that is in some sense similar to the one obtained 
in [8j. First of all let us recall Van der Corput's lemma(see for example [12], p. 332). 

Lemma 3.1. (Van der Corput) Let k > 1 be an integer, and <fi : [a, b] — > R such that 

\(j)( k \x)\ > 1 for all x G [a, b], and 4>' monotone in the case k = 1. 

Then 



b 



e 



it<t>{x) i)(x)dx 



b 



<c fc |t|-»(|HU- (0 , 6 )+ / |V(OIO, vt^o 



A first improvement has been obtained in [8] where the authors analyze the smoothing 
effect of some dispersive equations. We will present here a particular case of the results in 
[8], that will be sufficient for our purposes. In the sequel Q will be a bounded interval. We 
consider class A2 of real functions <fi G C 3 (Q) satisfying the following conditions: 

1) Set 5 = {£ G Q : <p" = 0} is finite, 

2) If £ ^ S<f, then there exist constants e, Ci, c 2 and a > 2 such that for all |£ — £ | < e ; 

Cl ie-eor" 2 <i0"(oi<c 2 ie-eor- 2 , 

3) 0" has a finite number of changes of monotonicity. 
Lemma 3.2. Let Q be a bounded interval, <fi G Ai and 

Jn 
Then for any x, t G R 

(3.1) \I{x,t)\<c^\t\- 1 ' 2 , 

where c^ depends only on the constants involved in the definition of class A2- 

Remark 1. The results of [8] are more general that the one we presented here allowing 
functions with vertical asymptotics, finite union of intervals or infinite domains. 

As a corollary we also have [8]: 

Corrolary 3.1. If <f> G A2 then 

e ^(O-^| " (O |i/ 2 ^ m <c^\t\- l i 2 (\\4,\\L~ m + /V(£M 

n v Jn 

holds for all x, t G R. 
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In the proof of our main result we will need a result similar to Lemma 13.21 but with 
|p///|i/3 i ns t eac i f Ip"] 1 / 2 in the definition of I(x,t). We define class A3 of real functions 
G C A (Q) satisfying the following conditions: 

1) Set S4 = {£ G Q : 0'" = 0} is finite, 

2) If £ G S^ then there exist constants e, c 1; c 2 and a > 3 such that for all |£ — £ | < e, 

(3.2) c^-U a - 3 < |0'" (01 <c 2 |£-£or 3 , 

3) 0"' has a finite number of changes of monotonicity. 
Lemma 3.3. Let Q be a bounded interval, G A3 and 

I(x,t)= [ e^^-^lf'iOl 113 ^. 



Then for any x, t G R 

(3.3) l/^^l^cjtr 1 / 3 , 

where c^ depends only on the constants involved in the definition of class A3. 

In the following we will write a < b if there exists a positive constant C such that a < Cb. 
Similar for a > b. Also we will write a ~ b if C*i& < a < C 2 6 for some positive constants 
C\ and C 2 . 

Proof. We observe that since f2 is bounded we only need to consider the case when t is 
large. 

Case 1: 0<m< \<f>'"{£)\ <M. 
We apply Van der Corput's Lemma with k = 3 to the phase function 0(£) — xt^/t and to 
V>= |0T /3 . Then 

\I(x,t)\ < C(tm)-i(||V|U»(n) + II^'IUh"))- 
Since 0"' has a finite number of changes of monotonicity we deduce that <p^ changes the 
sign finitely many times and then 

W\Ww = \l \(f(0rh {4) (0\d£ < \m-i I |0 (4) (OK < C(m,M). 

Hence 

\I{x,t)\ < C(M,m)ri 

Case 2: < |0'"(OI < M - 
Using the assumptions on we can assume that there exists only one point £0 G fi such 
that 0"'(£o) — 0. Notice that if G .A3, then any translation and any linear perturbation 
of (i.e. 0(£ — £0) + at; + b) is still in A3 and the conditions in the definition of set A3 are 
verified with the same constants as 0. Therefore we can assume that £0 = and 0'(£o) — 0. 
Moreover let us assume that as £ ~ 0, |0'(£)l ~ l£| a an d \4>'"(0\ ~ l£l f° r some numbers 
a > 2 and (3 > 0. 

We distinguish now two cases depending on the behavior of 0' near £ = 0. If a > 4 then 
|0( fc )(£)| ~ |£| a ~ fe as £ ~ for A; = 2, 3 and, in particular (3 = a — 3. The case a = 3 cannot 
appear since then /3 = a — 3 and does not vanish at £ = 0. For a = 2, |0'(£)l ~ |£|> 
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10" (01 ~ 1 as £ ~ and the third derivative satisfies |0'"(OI ~ l£l /3 as £ ~ for some 
positive integer /3. This last case occurs for example when 0'(O = £ + 0- m all cases 
(3>a-3. 

We split Q as follows 



I(x,t) 



e *W«)-^)|0'"(^)|5^ + / e^^-^10'"^)^^ = h + h- 

Since £ = is the only point where the third derivative vanishes we have that outside 
an interval that contains the origin <f)'" does not vanish. Thus 7 2 can be treated as in the 
first case. 

Let us now estimate the first term I\. We define fij, 1 < j '< 3, as follows 



Q, 



n x = {^ G H||e| < mm(e, Itr 1 /")}, 



£eft-fii||£|<e, and 0'(£) - 



X 


1 


X 





< - 





t 


- 2 


t 



fi 3 = {£e^-(^iUfi 2 )||£|<e}. 

In the case of Qi we use that for some j3 > 1, the third derivative of satisfies ci |^|^ < 
/"(OI<c 2 |el /3 for|e|<e. We get 






I0 w (oi^e <c| 

Hi Jf2i 

where the last inequality holds since a < /3 + 3 and \t\ > 1. 

In the case of the integral on f2 2 we assume that rr 7^ since otherwise f2 2 has measure 
zero. Observe that for £ e f2 2 we have 



±I0'(OIt 



< 



0<«) - f 



1 
< - 

- 2 



which implies that 



< 10' (01 < 5 



Since |0'(OI ~ IC| Q_1 we have that |f | ~ |x/t|^. Then |0'"(OI ~ ICP ~ k/t|^ an d 

'(01 >o. 



mm 

£en 2 



Applying Van der Corput's Lemma with k — 3 and using that 0' ' changes the sign finitely 
many times we obtain that 



9-2 



*<*<0-*y"(£)|*de < C(rnin|0'"(OIM)-^lll0'"(O|5||^(n 2 ) + IK^OI^'IU^ 



^en 



= C7(min|0'"(O|)-^|^(max|0'"(O|Ui/ |0 W (OI^ |0 (4) (OK 
<L7(min|0'"(O|)^max|0'"(O|^ri 
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Since on Q 2 , 10 (01 ~ \ x IA a ' 1 i there exists a positive constant C such that 

max|0'"(O|^<C(min|/'(O|)i 

which gives us the desired estimates on the integral on f2 2 . 

Now, we estimate the integral on fi 3 . Observe that we have to consider the case |t|~ 1//a < 
e, otherwise f2 2 = ^3 = 0. In particular, for £ G ^3, we have |t|~ 1//a < £ < e. Integrating 
by parts the integral on fl 3 satisfies 

(3.4) 

1 



e ^(e)-*O|0'"(£)|i^ 



n :i 



n,, 



>iMt)-xt)y ^'"(Ol 3 d t 



< 



±e 



-.wo-*) l^"(OI* 



0'(£)-flan 3 

^-^ li^oi-M^o^o - f)- 1^(01^(0 ^ 



0,3 



, 2 |<H0I* 

< -- max ■ 



1*1^3 u'(o-f mJn 3 w(0-7 



(0'(O-f) 2 

l<T(OHl0(4)(OI 



W"(0\^"(0\ d ^ 



1 

r ^-/fi 3 (0'(C)-f) 



In the following we obtain upper bounds for all terms in the right hand side of (j3.4p . 
Since on fl 3 , \<f/(£) — x/t\ > |x/2i|, there exists a positive constant c such that 



x 



4>'{t)-j >c\<p\o\>m a -\ v£Gfi 3 . 



In the case of the first term 
(3.5) 



1 mo 1 

sup 



\ C |£|3 

< — sup 



C 



sup 10 



l*U^ 3 |0'(O-f| ~ N^OaiO"- 1 l*Uen 
since |0 < e < 1 and |£p/3-a+i < |^|(«-3)/3-«+i = |£|-2«/3 < | t |2/3 



f-a+l < U.-1/3 



The second term satisfies 

1 r §i<r(or^ (4) (Ol 

n 3 



k 



WO-fl 






3 



|]SV 4) (0K<£ /_ IO^- Q+1 I0 (4) (OMO 



t 



Integrating by parts, applying the triangle inequality and using the definition of f2 3 we get 



10 3 -+v (4) (OI^<su P |0 



-f£-a+l| fjJUftW 1 / Itl^-'M ./■'■' 



n.s 



^3 



(01+ / \z\* a w"{zm 

o 3 



< 



sup 10- 



-a+l 



ft 3 



io* x 



fia 



< SUp |0 3 Q+1 < |t| 2 / 3 , 

n 3 
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where the last inequality follows as in (j3.5p . 

The last term in (13.41) can be estimated as follows 






Putting together the estimates for the terms in the right hand side of ( 13. 41) we obtain that 
the integral on Q 3 also decays as It]" 1 ^ 3 . 

The proof is now finished. □ 

4. Proof of the main result 

In this section we prove the main result of this paper. In order to do this, we will follow 
the ideas of [T] in the case of a discrete operator. Let us consider the system 

" iU t + AU = 0, 



(4,1; 1 U(0) = V , 

where U(t) = (u(t,j))j^ and operator A is given by f |2.Tj) . We compute explicitly the 
resolvent (A — A/) -1 , we obtain a limiting absorption principle and finally we prove the 
main result of this paper Theorem 11.11 

4.1. The resolvent. We start by localizing the spectrum of operator A and computing 
the resolvent R(X) = (A — A/) -1 . We use some classical results on difference equations. 

Theorem 4.1. For any b\ and b% positive the spectrum of operator A satisfies 

(4.2) a(A) = [-4max{^ 2 ,^ 2 },0]. 

Proof. Since A is self-adjoint we have that 

o~(A) C [ inf (Au, u), sup (Au, u)]. 

II U H;2(Z*)<1 IMIj2( Z »)<i 

Explicit computations show that 

(AU, U) = -b^ 2 ^ ( U j ~ U 3-l) 2 - tf + tf ( U -l ~ M 2 ~ b 2 2 ^2( U J+1 - U j) 2 - 
j<-l 1 2 j>l 

It is easy to see that (Au, u) < and 

(Au, u) > -2 maxj^ 2 , b 2 2 } ^ (u 2 + u 2 +l ) = -4 maxj^ 2 , b 2 2 } ^ u 2 . 

In order to prove that the spectrum is continuous we need to prove that for any A G 
[— 4max{65~ 2 , b^ 2 }, 0] we can find u n G Z 2 (Z*) with ||ii n ||i 2 (z*) < 1 such that ||(-A— \I)u n \\p(%*\ 
tends to zero. To fix the ideas let us assume that b 2 < 61 and A G [—Ab^ 2 , 0]. We construct 
u n such that all its components u n j, j < — 1, vanish. Thus for such w n 's we have that 

(Au n )j = b2 2 (A d u n )j, j > 1. 
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Using the fact that any A G [— 46^" 2 ,0] belongs to cr(6^" 1 A ( i) we can construct sequences 
( u n,j)j>i sucn that ||w n ||;2(z«) < 1 and \\(A — XI)u n \\p(2,*) — > 0. This implies that A G cr(A) 
and the proof is finished. □ 

Before computing the resolvent (A — XI)^ 1 we need some results for difference equations. 

Lemma 4.1. For any A G C \ [—4, 0] and g G I 2 (If), any solution f G Z 2 (Z*) of 

A d f(j)-Xf(j)=g(j), jVO 

with /(0) prescribed is given by 

(4.3) f(j) = tw-W + 2r _ 1 2 _ A £ r^*) 

feez* 

where a is determined by /(0) and r zs #ie unique solution with \r\ < 1 of 

r 2 — 2r + 1 = Ar. 
Moreover 

f(J) = Mr lj] + — ^ V(r^ fe l - r^+l fc l)a(A;), j ^ 0. 

k 

Proof. Let us consider the case when j > 1, the other case j < — 1 can be treated similarly. 
Writing the equation satisfied by / we obtain that 

f(j + 1) - (2 + X)f(j) + f(j - 1) = g(j), J > I- 

This is an inhomogeneous difference equation whose solutions are written as the sum 
between a particular solution and the general solution for the homogeneous difference 
equation 

f(j + 1) - (2 + X)f(j) + f(j - 1) = 0, j> 1. 
Let us denote by r\ and r 2 , |ri| < |r 2 |, the two solutions of the second order equation 

r 2 - (2 + A)r + l = 0. 

Since 2 + A G C \ [—2, 2] we have that r\ and r 2 belong to C \ R and more than that 
[r*i | < 1 < |r 2 |. Thus we obtain that 

(4.4) f(j) = ar{ + /^ + - -1 _ - £ r^'a(A;). 

fcez* 

Since / is an / 2 (Z + ) function we should have /3 = 0. Then formula (14. 3 p holds. The last 
identity is obtained by putting j = in (14. 4p and using that 2r — 2 — A = r — r _1 . □ 

As an application of the previous Lemma we have the following result. 

Lemma 4.2. Set Z 1 = Zn(-oo, -1] and Z 2 = Zn[l,oo). For any X G C\[-4max{o:f 2 , b^ 2 }, 0] 
and g G Z 2 (Z*) ; any solution f G / 2 (Z) of 

b: 2 A d f(j)-Xf(j)=g(j), jeZ s , 
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with /(0) prescribed is given by 

(4.5) f(j) = a ,rf + ^ ^\_ w E H^W*), J eZ s ,se {1,2} 

s s kez a 

where for s G {1, 2} ; constant a s is determined by /(0) and r s is the unique solution with 

\r s \ < 1 of 

r 2 s - 2r s + 1 = Xr s b 2 s . 

Moreover 

(4.6) f(j) = /(0)r?> + -^^(r^-rfl")^,] G Z.. 

r * r * fcez s 
The proof of this lemma consists in just applying Lemma [4. II to the difference equations 
in Z\ and Z 2 . 

Lemma 4.3. Let A G C \ [— 4max{o]~ 2 , b 2 2 }, 0]. For any o G Z 2 (Z*) t/iere exists a unique 
solution f G / 2 (Z*) o/ £/ie equation (A — XI) f = g. Moreover, it is given by the following 
formula 

_ b'l 
(4-7) /(j) = -^ ^- r [ E rftK*) + £ ri fc| a(fc) 

r s — r s x z — ' 
s s kez B 

where for s G {1,2}, r s = r s (X) is the unique solution with \r s \ < 1 of the equation 

r 2 -2r s + l = Xb 2 r s . 

Proof. Any solution of (A — XI) f = g satisfies 

A d f(j)-b 2 Xf(j)=b 2 g(j), jeZ s , 

bf(f(-l)-f(0)) = b 2 2 (f(0)-f(l)), 

where /(0) is artificially introduced in order to write the system in a convenient form that 
permits us to apply Lemma 14.21 
Using (14.61) we obtain 

/(-l) = /(0)n - b 2 £ r[ kl g(k) 

kez 2 
and 

f(l) = f(0) r2 -b 2 2 Y,r[ kl g(k). 

kez 2 
The coupling condition gives us that 

o x (l-nj + fej, (1 - r a ) s=1) tl 6Zs 
Introducing this formula in (14. 6p we obtain the explicit formula of the resolvent. □ 
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4.2. Limiting absorption principle. In this subsection we write a limiting absorption 
principle. From Lemma 14.31 we know that for any A G C \ [— 4max{&7~ 2 , b^ 2 }, 0] and 
ip G 1 2 {Z*) there exists R(\)(p = (A - A)~V £ ^( z ) an d it is given by 

(4.8) 

_ lil 

6 2 (l-r 2 ) + 6 x (l-ri)L^ ^ 

where r s = r s (X), s G {1, 2}, is the unique solution with \r s \ < 1 of the equation 

r 2 s -2r s + l = \b 2 r s . 

Let us now consider / = [— 4max{&7" 2 , b^ 2 }, 0]. As we proved in Theorem 14.11 we have 
that cr(A) = I. For any to G I and e > let us denote by rf e the unique solution with 
modulus less than one of 

r 2 — 2r + 1 = (u ± ie)b 2 r. 

Denoting r+ e = exp(;z+ e ) with zf e = af e + iaf e , af e < and af e G [— 7T, n] we obtain by 
taking the imaginary part in the equation satisfied by rf € that 

(exp(a+J - exp(-a+J) sin(a+ e ) = eb 2 s . 

Thus af e G [— 7r, 0]. A similar result holds for r~ t , a~ t G [0,7r]. 

Let us set r^ = lim e ^ rf e . Using the sign of the imaginary part of rf t we obtain that rf 
are the solutions with Im (rf) < < Im (r~) of the equation 

r 2 — 2r + 1 = w6 2 r. 

Also, using that r~ e = r+ e we obtain r~ = rf. 

For any u G J = I\ {-46^ 2 , -46^" 2 , 0} and ip G l l (Z*) let us set 

_( r ±\\j\ |- 

b 2 (l-r^ + b, (l- ri )L^ ^ 



+ r± _ ' t £((r±)*-*' - (^) lil+|fe| M^ 3 e Z s 



r — r 



We will prove that R (cj) are well defined as bounded operators from / 1 (Z*) to Z°°(Z*). 
We point out that we cannot define R ± (u) for u; G {— 46T" 2 , — 46^" 2 ,0} since for u = we 
have 7"i = r 2 = 1 and for w = 46~ 2 , s G {1, 2}, we have r s = — 1. We also emphasize that 
R~{uj)ip = R + {u)Tp. This is a consequence of the fact that for any uj G I, r~{u) = rf(u>). 
Formally, the above operator equals R(u ± ze) with e = 0. We point out that as operators 
on Z 2 (Z*), R(u) ± ze) are defined for any uj G I but only if e ^ 0. 
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Lemma 4.4. For any tp G / 1 (Z*) operator exp(itA) satisfies 

(4.9) 
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e itA tp = -L f e it0J [R + (u) - R-(uj)]ipduj. 

ZlTT J j 



Proof. To clarify the ideas behind the proof we divide it in several steps. 

Step 1. Let 1 1 be a bounded interval such that I C I\. There exists a constant 

1 1 1 



(4.10) 



C(u) 



Ym G L ^ 



MV2 \cobj + 41 1 / 2 |o;6l + 4|V2 
such that for all ui G I\ \ {—4b^ 2 , —Ab^ 2 , 0} the following inequality 

\(R(u ± ie)tp)(n)\ < C(w)|M|ji(z*), for all tp G l\Z*) and n G Z*, 

holds uniformly on small enough e. 

Step 2. For any we J, i? ± (a;) are bounded operators from / 1 (Z*) to Z°°(Z*) and 

ll-R (w)||ji(z*)-j°°(z*) < C(w). 

Step 3. For any u G J, ip G / X (Z*) and n G Z* the following holds 

lim(i?(w±ie)^)(n) = (if^w^Xn). 

Step 4. For any tp G / : (Z*) and n G Z* we have 

lim [ e ituJ (R(cu±ie)tp)(n)dcu= I e it0J (R ± (u)tp)(n)duj. 



Step 5. For any </? G /*(Z 



i/"y*-\ 



e iM p 



1 
2i?r 



e ltw [R + (u)-R-(uj)]tpduj. 



Proof of Step 1. Observe that for any wGl and e > we have 

\(R(u±ie)<p)(n)\ 

1 1 



< 



\ b - 2 (l-rt) + bT 2 (l-rt 



Solution rf t of equation r 2 — 2r + 1 = (a; ± ie)b 2 s r satisfies 



r ie-( r ty i \ 



i r v _ ( r i) _i i 



,.± 



-\r ± \ < 



r* - — 



6Jw±ie| 1/2 . 



Then for all u E I± and e small enough we have 

2 



|r ± I > 

I s.el — 



b 9 \u±ie\ 1 / 2 + (V i s \u>±ie\+4) 1 / r < 



> C>0 



and 



|r ± I < 



< C 1 < oo. 
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Thus for any ui G Ij 


we have 














\ r ± 

1' s,e 


1 

- try- 1 ! ~ 


1 

|l-r±J|l + r±J 


< 


|1 


1 


7 + 


1 

|l+r±| 


Using the equation 


satisfied by rf 


\ we find that 














|l-r±|=6 s 


|o;d=ze| 1 /2[ r ± | > 


\uj 


± 


iel 1 / 2 


> | 


col 1 / 2 



and 

|1 + r ± £ | = |(w ± ie)b 2 s + 4| 1 / 2 |r s ± J > |(w ± ie)b 2 s + 4\ 1/2 > \ub 2 s + 4) 1 / 2 . 

Putting together the above estimates for the roots rf e we find that for all uj E I\ and e 
small enough the following holds 

1 111 1 

~i T / -I- n 77 ~r 1 T / -I- \ 77 ^, "i iTTo ~r "i To : TTTTn "T 



i',., w L( , , ^-(rJJ-^^H 1 / 8 M? + 4|V2 | w6 | + 4H/2- 
We now prove that 

1 <^_ 



„± ■ 

2,e. 



+ ^T a (l-r£ e )|~ M 



We recall that the sign of the imaginary parts of rf e and rf e is the same. Also, since 
\rf e \ < 1, the real parts of 1 — rf e and 1 — rf e are positive. These properties of the roots 
imply that 

I& 2 - 2 (1 - r%) + bf(l - rfj\ > bf\l - rfj + bf\l - r± | > \uj\ 1 ' 2 . 

Putting together the above results we obtain that Step 1 is satisfied with C{uj) given by 
(I4.10p uniformly on all e > sufficiently small. 

Step 2 follows as Step 1 by putting e = and replacing rf e with rf. 

Proof of Step 3. We write 

R{ui ± ie)(p(n) = N^ R(u ± ie, n, k)cp(k), 

fcez* 

where R(u ± ie, n, k) collects all the coefficients in front of <p(k) in formula (14. 7\i . 

Using that, for any uj G J, rf e (u>) — > rf(u>) we obtain that R(u> ± ie,n,k)ip(k) — > 
R ± (u),n, k)(p(k). Since for any u G J and e small enough we have the uniform bound 

\R(u ± ie,n, k)(p(k)\ < C(u)\(p(k)\,Vk G Z*, 

we can apply Lebesgue's dominated convergence theorem to conclude that 

y R(u ± ie, n, k)<p(k) — > 2_, R (w, n, k)cp(k), 

fcez* kez* 

which proves Step 3. 

Step 4 follows by Lebesgue's dominated convergence theorem since we have the pointwise 
convergence in Step 3 and the uniform bound in Step 1. 
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Proof of Step 5. Applying Cauchy's formula we obtain that 



ua = J_f e itw R(u)dw 
2Wr 



for any curve T that rounds the spectrum of operator A. For small parameter e we choose 
in the above formula path T e to be the following rectangle 

T e ={tu ± ie, to G [— 4max{6^~ 2 , b^ 2 } — e, e]} 

U {— 4max{6]" 2 ,6^ 2 } — e + ir],r) G [— e, e]} U {e + ir],r] G [— e, e]}. 

Using the estimates for R(X), A G r e obtained in Step 1 and the convergence in Step 4 we 
obtain that for any if G / X (Z*) the following holds: 

e «A = J_ f e ittJ (R + (u;) - R-(u))<pdu. 
2m J j 

The proof is now complete. □ 

4.3. Proof of the main result. We now prove the main result of this paper. 
Proof of Theorem ! 1.1[ For any ip G / X (Z*) Lemma [4.41 gives us that 

(e itA <p)(n) = — I e itu (R + (u) - R-(u))cp{n)ds, n G Z*, 
2m J j 



where / = [— 4max{6 1 ,b 2 },0]. Using the fact that R (u)ip = R + (lo)<p we obtain 

(e ltA ip)(n) = - [e ituJ ((ImR + )(tu)cp)(n)du, n G Z*, 
71 J i 

where Im R + is given by 

(R+(u)<p)(j)-(R-(u)<p)(j) 



Im R + )(u>Mj) 



2i 
___ -( r +)\j\( r +)\k\ 



^"(l-rj)+^(l-rf) 

+ E ^(fe)Im & ; 1 ((r+)l^ fc l - ( r +)WI+l*l), j G Z s 

and for s G {1,2}, r+ is the root of r 2 — 2r + 1 = ub 2 s r with the imaginary part nonpositive. 

In order to prove fll. 12j) it is sufficient to show the existence of a constant C = C(bi, b 2 ) 
such that 
(4.11) 

^^ f ( r +)\J\(r+)\ k \ 

£ 1^)1 1 //^ ^(i -r+) + 6r 2 (l ~ ri) dU \ ~ ° m + 1] IM| W' Vj G Z *' 
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(4.12) £>(*)! I fe^lm 

k£Z s ^ I 



( r +)\j-k\ 

-duj 



r 7 — r: 



-i 



<C(|i| + l)-^|Mhz*), VjeZ*. 



The estimates for the other two terms occurring in the representation of Im.R + (u) are 
similar. 



Step I. Proof of (Q2]) . We prove that 



(4.13) 



sup 
j'ez 



ilil 



e ltw Im ■ 



-gL> 



■r*~r I 'T ,_ l 



<C(6 1 ,6 2 )(|«| + 1)- 1 ^ 3 , VtGM. 



We split I as J = ii U I 2 where J x = [-4 max{frf 2 , b 2 2 }, 4bj 2 } and I 2 = [4&<T 2 , 0]. If u G J u 
the following equation 

1 



has real roots and then 



r + - = 2 + ubi 
r 

( r +)\j\ 
e ^ i m v >> du = 0. 



When uj G J 2 , root r s of equation r s + — = 2 + w& 2 has the form r s = e l0 , 9 G [0, n] 
Using the change of variables u = b; 2 {2 cos 8 — 2) we get 



\\j\ 



e ttuJ Im ■ 



- fr+W 



du = 2K 2 / e ^" 2 (2cos e - 2 ) Im . 



-W 



o 



e lK > - e l 



sinOde 



(,-i\j\e 
2b; 2 I ^^^-^Im^—^siufhif) 



o 



2z sin 6 



__ b -2 / e ^ 2 (2cose- 2 ) Ree - i |i|^ 

e itbr (2co86,_2) (e i|j ' |fl + e~ i|3 ' |fl )d0. 



2 Jo 

Van der Corput's Lemma applied to the phase function 

that 



(2 cos 9 - 2)bZ 2 +j9/t shows 



(4.14) 



e it(2cose-2)b s 2 e ij6 d Q 



< C(b s )(\t\ + I) -3 , V t G M, Vj 6 Z 



The proof of (I4.12p is now finished. 

Step II. Proof of (14.111) . It is sufficient to prove that 



sup 



„itw 



[rT ) J (U 



6iT a (l-rJ) + ^ 2 (l-r+) 



do; 



<C(6 1 ,6 2 )(|t| + l)- 1 / 3 , Wg 



To fix the ideas let us assume that 6 2 < &i • We split interval / as follows I = 1% U J 2 where 
J x = [_4fe- 2 , -46]" 2 ] and J 2 = [-46^ 2 , 0]. We remark that on h, rf G R and r+ G C \ R. 
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On I 2 both rf and r 2 belong to C \ R. We prove that 



(4.15) 

and 
(4.16) 



sup 

j,fceN 



Atui 



(rt) j (r+) k 



h b 2 \l-rt) + b^(l-rt 



-duj 



<C(6 1 ,6 2 )(|t| + l)- 1 / 3 



sup 

i,fceN 



AtiU 



(rt) J (4) k 



-du 



<c{b 1: b 2 ){\t\ + ^)- xlz - 



h 6 2 - 2 (l-r+)H-for 2 (l-^) 

Let us set h(u>) = b 2 2 (l — r 2 (u)) + b~ 2 (l — rf(u>)) Using the same arguments as in the 
proof of Lemma S3 we get that \h(u)\ > C(6i, b 2 )\co\ 1/2 . Then, on h, \h{u)\ > c > 0. 
Moreover |/i'(u;)| — c 2 < °°- Using integration by parts we obtain that 



Stiii 



(rt) j (r 2 



h b?(l-r+) + b?(l-rt 



-du 



< sup 



„ito / 



I/;., 



rfy(rt) k du(\\l/h\\ L ~ {h) + \\(l/h)% Hli: 



< C(bi,b 2 ) sup 
A similar argument shows that 



„itw 



{r+y{rt) k duj 



e itu> (rf) j (r+) k du 



16" 



< sup 



16.; 



46- 



e^(r 2 + ) fc ^(||(r^||^ (7l) + ||((r+y)'|| L . (/l) 



Observe that for u E I±, r^(u) given by 

+ . 2 + b\u-^{2 + b\uY-A 
r i» = g 

is a decreasing function. Thus 

\\((rtyy\\L Hh ) < || WlU^) < 1, VjGN. 
The proof of (I4.15P is now reduced to the following estimate: 



sup 
3/e/i 



e ita; (r 2 (w 



)) fc cia; < C(6i,6 2 )(|t| + 1)- 1/3 ,VA; 6 N,t e 



16.T 



Making the change of variables to = b 2 2 (2 cos 9 — 2) and applying Van der Corput's Lemma 
as in the final step of Step I we obtain that 



e ltui (r+(u)) k du 



167 



2bn 



Jt%(2co S 6-2) e -ik6 sin9du 



2arcsin(62/?/) 



< C(6 2 )(|t| + 1)-V3_ 



We now prove ( 14. 16ft . We first make the change of variables u> = b 1 2 (2 cos 9 — 2). Thus 



Atui 



(rt) j (4) k 



du = 2b~ 2 I c itb ~ 2 (2 cos e ~ 2) c~ ije c~ 2ik axcsin ( b ^ b i 1 sin f ) 01 " u d9 



sin 9 

W) d0 



/ 2 6 2 - 2 (l-r 2 + ) + 6r 2 (l-r+) 
where h(9) = b 2 2 (l -r£(6)) +6^ 2 (1 -rf (0)), rf(9) = e~ ie andr+(#) = e-MaxcrinCWf 1 ■»*!). 
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Using that far from 6 = function h satisfies \h(6)\ > we choose a small parameter e 
and split our integral as follows: 



itbT 2 (2cose-2) -ijO ~2ika,Tcs'm(b 2 b~ 1 sin 4) S1D r lQ _ T -I- T 

h(0) 



,„^7 2 (2 cos 6-2) -ij9 -2ik arcsinfefer 1 sin §) Sm " ^fl 



„«<&, 



jSin^ 

w 



sinO 



, / Jib' 2 (2 cos 6-2) -ij6-2ik arcsinfefc," 1 sin f) bm ^ r/i 

- / e ^at/ 



Observe that on interval [0, e] 

sin# 



h(6) 



L°°(0,e) 



+ 



d ^in^, 



and on interval [e, it] 



HO) 



L°°(e,7r) 



+ 



d9 y h(9) 



d , 1 



Li(0,6) 



< M < oo 



de y h(e)' 



L x (e,-K) 



< M < OO. 



Then we have the following estimates for 7\ and T 2 



and 



I Ti| < M sup 

xe[0,e] 



|T 2 | < M sup 

x£[e,7r] 



^iib-L (2cos0— 2) — ij'S — 2ifcarcsin(b2fr 1 s ^ n 2^rlft 



e itb 1 2 (2 cos 8-2) e -ij0 e -2ika,rcsm(b 2 b 1 1 sin|) g j n ^^ 



We now apply the following lemma that we prove later. 

Lemma 4.5. Le£ a G (0, 1] and < 5 < it. There exists C(a,5) such that for all real 
numbers y, z and t 



(4.17) 

and if S > 
(4.18) 
We obtain that 

and 



e ii(2 cos 0+22 arcsin(a sin §)) e ij/0 ^ q^q 



< C(a,<5)(|t| + I)" 1 / 3 



it(2 cos 8+2z arcsin(a sin ^)) iyd in 



< C( a ,5)(|t| + 1)- 1 / 3 . 



|T 1 |<MC(a,e)(|t| + l)- 1 / 3 



|T 2 |<MC(a,e)(|t| + l)- 1 / 3 . 
The proof of Theorem 11.11 is now finished. 



a 
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Proof of Lemma \4-5\ Since the integrals in f 14 . 1 Tj) and (14.181) are on bounded intervals it is 
sufficient to prove that, for t large enough, each of the integrals is bounded by |t| -1 ^ 3 - In 
the case of (14.171) we will consider the case 5 = since the proof for 5 > is similar. 
Let us denote by ip either the function X(o,it-S) or si n $• We set 

n 

p(9) = 2cos# + 2zarcsin(asin -), 9 6 [0, ir]. 
Using the Maple software we obtain that 



mm [6W + WW > ™ {4 + ^=^-, ^{, ~ ^f}. 



If z is such that \z — a | > e > then Van der Corput's lemma applied to the phase 

function p{9) + y6/t guarantees that 



e itm e iye i>{9)d9 < C{a,e)(\t\ + 1)" 1/3 . 
o 



Assume now that \z — A ^\ a \ < e with e small enough that we will specify later. Let us 
write 

4v / r 3 ^ 

z= Vb 

a 

with b a small parameter such that \b\ < e. With this notation p(9) = pb{9) = q(9) + br(9) 
where 



,/a ^VT^aJ . 9^ 

q(9) = 2 cos{9) -\ arcsin(a sin -) 

and 

o 
r(9) = 2arcsin(asin -). 

Solving system (q"(9),q'"(9)) = (0,0) with Maple software we obtain that it has a unique 
solution 9 = tv. Thus for any 5 < n there exists a positive constant c(a, S) such that 

\q"{9)\ + \q'"{9)\ > c{a,5), \/9e[0,tt- 5]. 

It implies the existence of an e = e(a, S) such that for all |6| < e 

\d(d)\ + \p'm\ > <a,5) - \b\ sup (\r"\ + \r'"\) > M, V 9 E [0, tt - S). 

xe[0,7r] * 

Hence, Van der Corput's Lemma applied to the phase function Pb{9) + y6/t guarantees 
that 

/ e itpb{e) e iye ^{9)d9 < C{a,8)(\t\ + 1)~ 1/3 , V|6| < e,V t,y E E. 
Jo 

The proof of (I4.18P is finished. 
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To prove estimate (I4.17P it remains to show that we can choose 5(a) small enough such 
that for all \b\ < e 

(4.19) \h(t)\:= / e itpb{e) e iye sm(6)d6 <C(a)(\t\ + l)- 1/3 , Vy,teR. 

The Taylor expansions of q and r near 9 = tx are as follows 

fn . -2a + 8v / T^^arcsin(a) 1 (2 a 2 - 1) (9 - vr) 4 1 (4 a 2 - 1) (9 - vr) 6 

gW = a 16 ^TT^ 384 (_l + a 2 ) 2 + ° (( ^ 7r) 

and 



v2 , 1 a (2a 2 + l) 4 



r(0) = 2 arcsin (a) - - -=== (0 - tx) 2 + — - v ,1/ (0 - vr) 4 + O ((9 - vr) 6 ) . 

y J ^ ' 4 v^^^2 ^ ^ ^ 192 (1 _ a 2)3/2 ^ ' T W M 

Also the second derivatives of g and r satisfy 

3 (2a 2 -l)(g-7r) 2 4 

g(0) = -- _ 1 + a 2 + 0(|fl-vr| ) as^~7r, 

and 

r"(^) = -i-=^ = = + 0(^-7r) 2 as^vr. 
^ V 1 — o 

Observe that for a ^ \j\f2, the second derivative of g behaves as (0 — tx) 2 near 9 = tx. 
Otherwise it behaves as (9 — tx) a near the same point. Since the proof of (I4.19P is quite 
different in the two cases we will treat then separately. 

In the sequel 5(a) is chosen such that we can compare q and r with their Taylor expres- 
sions near 9 = tx. 

Case 1. a^ l/v2. The main idea is to split the interval [tt — 5(a), it] in three intervals 
where we can compare 1 6* — -zr | with jo) 1 ' 2 and decide which of them dominates the other: 

[vr - 5(a), vr] = [vr - 5(a), tx - a 2 \b\ 1/2 } U [vr - a 2 \b\ 1/2 , tx - «i|6| 1/2 ] U [vr - a x \b\ l/2 , tx], 

where a± « 1 << a 2 are independent of b but depend on the parameter a. More precisely 
the parameters a± and a 2 are chosen in terms of the first two coefficients of the Taylor 
expansion of functions q and r near 9 = tx. 

Let us consider the interval [it — 5(a), tx — a^lo) 1 / 2 ] with a 2 large enough. In this interval 
1 6* — tt | dominates l&l 1 / 2 and we apply Lemma [3721 We check the hypotheses of this lemma. 
In this interval the first derivative of p b is of the same order as 1 9 — tx | 3 : 

\p' b (9)\ > W(0)\ - \b\\r'(9)\ > Cl \9 - tx\(\9 - tx\ 2 - C 2 \b\) > C 3 \9 - vr| 3 

and 

\p' b (0)\ < W(0)\ + \b\\r'(0)\ > Ci\9 - tx\(\9 - n\ 2 + C 5 \b\) > C 6 \9 - tx\ 3 . 
Also, the second derivative satisfies: 

\p'l{0)\ > W'(0)\ - \b\\r"(9)\ > C 7 (\9 -n\ 2 -C 8 \b\)> C 9 \9 - tx\ 2 

and 

\p'l(9)\ < \q"(e)\ + \b\\r"(9)\ > C 10 (\9 - tx\ 2 + Cu|6|) > C l2 \9 - tx\ 2 . 
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We emphasize that all the above constants are independent of b. Observe that on the 
considered interval \p%\ > \b\. If we try to apply Van der Corput's Lemma with k = 2 we 
obtain 

»7r-a 2 H 1/2 



e ltpbW e iyy sm(9)d9 < (\tb\) 



-1/2 



ir— 5(a) 



max \ sin 6\<C (5 (a))\tb\- 1/2 , 

■k— 5(a), it— ct2 |f'| 1/ ' 2 ] 



an estimate that is not uniform in the parameter b. 

However, using Lemma T3.2I we obtain the existence of a constant C depending on all the 
constants Ci,i = 1, ..., 12 but independent of the parameter b, such that 

(4.20) 



7T — «2| 



11/2 



it— 5(a) 



e ltPb{e) e iye sm(9)d9 

^. . 1 /9 / I sin(t/J 

<C\t\- l ' 2 ( max ±— \r{ 



7T-Q 2 |f)| 1/2 



itp b (9) iydi n "(Q\\l/2 



>e iy0 \ P : 



sin(0) 



it— 5(a) 



/2 



IT — Q 2 | 



,|l/2 



it— 5(a) 



< c\t 



-1/2 



sin 



max 



< C\t\- 1/2 max 



\P'1(0)\ 1/2 

MO) \' (f} 
|#(*)|W l 

I sin 



dO 



dO 



<cir 1/2 . 



[■K-5(a),ir-a 2 \b\ 1 / 2 ] IK'WI^ 2 ' ' [7r-5(a),7r-a 2 |b| 1 / 2 ] \9 ~ 7r| 

On the interval [7r — a2|o| 1//2 ,7r — ai|o| 1//2 ] the third derivative of pb satisfies: 

|/'(0)|~|0-7r||C(a)+&|~|&| 1/2 , 

since C(a) 7^ in the case a 7^ l/\/2. Applying Van der Corput's Lemma with A; = 3 we 
get 

(4.21) 



7T — CK1 



|6|V2 



7T— a 2 



|6]V2 



i it ^^e iv9 wa{e)de < {\tb\ 1/2 y 1/3 



max 



in 01 < itr 1/3 . 



sin 



6e[n-a 2 \b\ 1 / 2 ,ir-a 1 \b\ 1 / 2 ] 



On interval \n — ai|o| 1,/2 , 7r] with cti small enough, the term \br"{9)\ dominates \q"{9)\. 
The the behavior of p'l{9) is given by \br"{6)\: 

\p'l{9)\ > \br"{9)\ - \q"{9)\ > Cl (\b\ ~ C 2 \0 - n\ 2 ) > C 3 \b\, 

for some positive constants C\ and C2 independent of the parameter b. Applying Van der 
Corput's Lemma with k = 2 we get 



(4.22) 



Tr-ailbp/a 



e itpb{e) e iye sm{9)d9 < (\tb\) 



-1/2 



max 

eel-ir-a^b] 1 / 2 ^] 



sm 



01 < itr 1 / 2 . 



Using PD| . 032U) and (Q2|l we obtain that (Q5]l holds uniformly for all \b\ < e,y 
and t real numbers. 

Case 2. a = 1/a/2. In this case the Taylor expansion of function q at 9 = tx is given by 

-2a + 8 v / T^^ 2 arcsin(a) 1 (4 a 2 - 1) [9 - vr) 6 



#) 



384 



-1 + a 2 



O(|0-7r| 8 ). 
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We split the interval [it — 5(a), it] as follows: 

[tt - 5(a), ix] = [ir - 5(a), vr - a 3 |6| 1/4 ] U [tt - a 3 |6| 1/4 , vr - a 2 |6| 1/4 ] 
U [tt - a 2 \b\ l/ \ vr - ai \b\ 1/2 } U [tt - a x \b\ 1/2 , vr], 

where a 2 << 1 << a 3 and all a±, a 2 , a 3 are independent of b. 

On the first interval [ir — 5(a), n — a^b] 1 ' 4 ] we apply Lemma T3.31 We have to check that 
the first third derivatives behave as powers of \9 — tt | in this interval. Observe that 



\p' b (o)\ > c x \e - n\(\e - tt| 4 - c 2 |6|) > c 3 \ 



71 



and 



\p' b (9)\ < Ct\9 - tt|(|0 - tt| 4 + C 6 |6|) > C 6 \9 - tt| 5 . 
In a similar manner 

C 7 |^-7r| 4 <|KWI<C 8 |^-7r| 4 . 

Also the third derivative satisfies 

\p'b(0)\ > C 9 \9 - tt\(\9 - tt\ 2 - C 1Q \b\) > C u \9 - vr| 3 

and 

\Vb(0)\ < C l2 \6 - tt\(\9 - vr| 2 + C 13 |6|) > C 14 \9 - n\ 3 . 

We now apply Lemma [3731 taking into account that all the above constants are independent 
of b and we obtain 



(4.23) 



7r-a 3 H 1/4 



e itPb{0) e iye sin 9d9 



Tr—S(a) 



TT-a 3 \b\ 1 / i 



itPb(d)JyO\„"/(a\\i/3 



>e^\p 



sin# 



tt— 5(a) 



<\t\- 1/3 ( m, 

V \tt— 5(a), 7r 
< \t\~ 1/3 



sin#| 



max 



TT— 03\b\ 



dhl 1 /" 



[rr-S(a)^-a 3 \b\^} Wb( e )\ 1/3 J*-S(a) 

I sin 9 1 



\pm\ i/3 

sin 9 

WW 3 



d9 



d9 



max 



rc-5(a),n-a 3 \b\V±] \p' b " (9)\ 1/3 



< i*r 1/3 



sin^l 



max 



TT-S(a),TT-a 3 \b\ i -/ i } \9 — 7r| 



< cm- 1 ' 3 . 



In the case of the interval [tt — a 3 |6| 1//4 , tt — a 2 |5| 1//4 ] we apply Van der Corput's Lemma 
with k = 3 and use that 



Then 
(4.24) 



K'(9)\>C 1 \9-n\ 

7T-a 2 \b\^ 4 



TT 



e itpb{e) e iy8 sin9 



TT-a 3 \b\ 1 / A 



C 2 \b\) >C,\9- n\(a 2 2 \b\ 1/2 ~ C 2 \b\) > CM 1 '^ 1 ' 2 - 



< (\t\\b\ 3/4 )- 1/3 max |sin^| <C|t|- 1/3 . 

[7r-a 3 |b| 1 /4 i7r _ Q , 2 | |i/4] 
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Let us now consider the integral on the interval \n — a2|&| 1//4 , vr — atil&l 1 ' 2 ]. Observe that 
in this case 



(4.25) 



w-axW 1 / 2 



e itpbW e iye sin 9d9 



TT — Ct2\ 



.11/4 



7T — Ctl| 



,|l/2 



< / \sin9\d9< / \sin6\d6 



a 2 \b\ 



1/4 



/■a 2 |b| 1/4 

< / ode < c\b\ 1/2 < c\t\ 



-1/3 



as long as \b\ < |£|~ 2 / 3 . 

We now consider the case |6| > |£| -2 / 3 and prove that a similar estimate can be obtained. 
Observe that on the considered interval the second derivative of pb satisfies 

Wl(0)\ > \b\\r"(9)\ - \q"(9)\ > d(|6| - C 2 \9 - tt| 4 ) > C x (|6| - C 2 (a 2 \b\ 1 ^) 4 ) > C 3 \b\. 

Thus, Van der Corput's Lemma with k = 2 gives us 

(4.26) 



Tr-ail&l 1 ^ 



TT — 02 



|6|V4 



e itpb e iye sin 9d9 < (\tb\) 
< \t\- 1/2 \b\-^ < \t 



-1/2 



max 



0e[7r-a 2 H 1/4 ,i— aiH 1/2 l 



sin#| < (Itfel)- 1 / 2 ^! 



-V2IM1/4 



1/2 n I 1/6 



Itl 



-1/3 



On the last interval [7r — «i l&l 1 / 2 , it] the term \br"(9)\ dominates \q"(9)\. Then the behavior 
of p'l{9) in the considered interval is given by \br"{ 



\p , b '(9)\>\br"(9)\-\q"(9)\>C 1 



Co\9 



71" 



>C 3 |6|. 



Thus 

(4.27) 



»r-ai|b|Va 



e itpb{e) e iy6 sin(9)d9 < (\tb\) 



1-1/2 



max 



sm 



0i < itr 1/2 . 



eeiTr-aiifip/z^] 



Using the previous estimates (Qgjl . f T4T2"4j) . (1Q5) . (j£25]) and (T4T2TJ) we obtain that 
estimate (14. 19j) also holds in the case a = l/\/2. 

The proof of Lemma 14.51 is now finished. □ 



In the case of system (12.8B the proof of Theorem 12.21 follows the lines of the proof of 
Theorem 11.11 by taking into account the representation formula for the resolvent of the 
operator A given by ( I2.9J) . 

Lemma 4.6. Let X G C \ [-4 max{6^ 2 , b 2 2 }, 0] and A given by (|2T9|) . For any g G / 2 (Z*) 
t/iere exists a unique solution f G Z 2 (Z*) o/ i/ie equation (A — XI) f = g. Moreover, it is 
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given by the following formula 

(4.28) 



m 



-r 



\j\ 



K (rT -n) + b2 2 (rj 1 -^) 



9(0) + E r i l 9(k) + Y, r 2 l 9(k) 



keZ! 



+ 



r s -r s 



n£< 



kez 2 



\j-k\ _ Jj\ + W 



r)>W\)g(k), jeZ s 



keZs 



where for s G {1, 2}, r s = r s (X) is the unique solution with \r s \ < 1 of the equation 

r 2 s -2r s + l = Xb 2 s r s . 
We leave the complete details of the proof of Theorem 12.21 to the reader. 

5. Open problems 

In this article we have analyzed the dispersive properties of the solutions of a system 
consisting in coupling two discrete Schrodinger equations. However we do not cover the case 
when more discrete equations are coupled. The main difficulty is to write in an accurate 
and clean way the resolvent of the linear operator occurring in the system. Once this case 
will be understood then we can treat discrete Schodinger equations on trees similar to 
those considered in [4] in the continuous case. 

There is another question which arises from this paper. Suppose that we have a system 
iUt + All = with an initial datum at t — 0, where A is an symmetric operator with 
a finite number of diagonals not identically vanishing. Under which assumptions on the 
operator A does solution U decay and how we can characterize the decay property in terms 
of the properties of A? When A is a diagonal operator we can use Fourier's analysis tools 
but in the case of a non-diagonal operator this is not useful. 
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